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A VECTOR STUDY OF LINEARIZED SUPERSONIC FLOW 
APPLICATIONS TO NONPLANAR PROBLEMS 1 

By John C. Martin 


SUMMARY 

A vector study oj the partial-differential equation oj steady 
linearized supersonic flow is presented. General expressions , 
which relate the velocity potential in the stream to the conditions 
on the disturbing surfaces , are derived. In connection with 
these general expressions the concept oj the finite part of an 
integral is discussed. 

A discussion of problems dealing with planar bodies is given 
arid the conditions for the solution to be unique are investigated. 

Problems concerning nonplanar systems are investigated , and 
methods are derived for the solution of some simple nonplanar 
bodies. The surface pressure distribution and the damping in 
roll are found for rolling tails consisting oj jour, six, and eight 
rectangular jins for the Mach number range where the region of 
interference between adjacent jins does not affect the jin tips. 

INTRODUCTION 

In the presentation of the theory of the flow of an idealized 
incompressible fluid, vector methods can be used to reduce 
greatly the mathematical manipulations involved. The 
study of steady linearized supersonic flow may also be aided 
by the use of vector methods. Two types of approaches, 
however, can be used. Perhaps the more obvious is to make 
use of common vector methods as was done in reference 1. 
The other vector method, which was introduced by Robinson 
in reference 2 and is used in this report, appears to be more 
suited to the study of the linearized partial-differential equa- 
tion of steady supersonic flow. This method allows a deriva- 
tion of a hyperbolic scalar potential and a hyperbolic vector 
potential along lines analogous to the derivation sometimes 
used (ref. 3, ch. VIII) in dealing with common scalar and 
vector potentials. 

The present report presents a vector derivation of many 
general results which have been found by various methods 
and are given in the published literature on the linearized 
partial-differential equation of supersonic flow and also 
presents some results which are not found in the literature. 
The general results of Hadamard (ref. 4, p. 207), Puckett 
(ref. 5), and Heaslet and Lomax (ref. 6) are found as special 
cases of a general expression for a scalar potential, and the 
results found by Robinson (ref. 2) are obtained by the use 
of a vector potential. The derivation of the scalar potential 
doubtlessly helps to clarify the concept of the finite part of 
an integral. 


A discussion of problems dealing with planar bodies im- 
mersed in a supersonic flow is given, and the conditions 
necessary for the solution to be unique are investigated. 

Problems dealing with nonplanar systems are also dis- 
cussed, and methods are derived for the solution of some 
simple problems dealing with nonplanar bodies. The surface 
pressure distribution, the spanwise loading, and the damping 
in roll are found for rolling tails consisting of four, six, and 
eight rectangular fins for the Mach number range where the 
region of interference between adjacent fins does not affect 
the fin tips. 


SYMBOLS 


A 

A 

a 

C\, Ci, Cz 
c 

E , F 

f 

G 

H 

6/2 

h j k 

M 

n 


hyperbolic vector potential 
aspect ratio of tail fin 
positive constant 
arbitrary constants 
chord 

arbitrary vector functions 
scalar function defined bv equation (19) 
vector function associated with vector 
function F 
vorticity vector 
span of tail fin 

unit vectors in x-, ?/-, and 2-directions, 
respectively 
Mach number 

unit vector normal to element of area da 


n,,= — i& 1 v\-\-jv2 J r kv z 
n h ' = — ipW +jW + kvj 
n h *=~i0W+jv 2 *+kv 3 * 

A C p pressure-difference coefficient 

p rate of roll 

Q function used in equation of surface of 

discontinuity 

q part of velocity vector which is made up of 

hyperbolic curl of vector potential 

q' total perturbation velocity 


R=J(x-i-) 2 -(S 2 (y-vY-l3 2 (2-0 2 

IF small constant 


r=^(x-zy+(ny-r,y+(Hz-n z 

S area of tail fin 

S 0 surface of discontinuity 


Supersedes NACA TN T 2641, “A Vector Study of Linearized Supersonic Flow Applications to Nonplanar Problems” by John C. Martin, 1952. 
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Si, S2, S3, 64, S3, 
S', T, T’ 

V 

r„, v u V 2 

J, ]h 2 


surfaces of integration 

free-stream velocity 
volumes of integration 


Cartesian coordinates (j-axis parallel to 
free-stream direction) 


p=^jAP— 1 

T spanwise circulation 

/ r t . e . xr r t . e . \ 


small positive quantity 


f, n, r 

Cartesian coordinates (£-axis parallel to 
f ree-st ream direction) 

8, p 


polar coordinates 

A , \p 


scalar functions 

A 


given volume 

Vu r 2 , 

V3 

direction cosines of outward normal to 
element of area da 

V\ , Vo 

\ v* 

direction cosines of normal (directed away 
from point (x,y,z)) to surface S 0 

Vl*, Vo 

*, *,* 

direction cosines of normal to element of 

area da used in equat ion (45) 

a 


slope of deflected area 

T 


area of integration 

<f>, <f>», </>!> 4>2> ^0 

scalar potentials 

Cl 


rolling-moment coefficient per fin, 
Rolling moment per fin 

1 v^ b 

2 p 6 2 


" »a, 1 

I 


af . ( fc/2)J 



Ip *0 

<£ 


indicates integration over closed line or 

j 


surface 


denotes finite part of integral 


THEORY 


This report deals with the linearized partial-differential 
equation of steady supersonic flow. This equation is given 
by 

— /3 2 <^xr+0//i/ + < / ) 2.’— 0 (0 


The potential is assumed to be continuous in the stream 
direction, and the potential is assumed to be always finite. 
Assuming the potential to be finite and continuous in the 
stream direction has the effect of requiring the aerodynamic 
lift and moment (calculated by use of the linearized pressure) 
of finite bodies to be finite since the linearized pressure is 
related to the derivative of the potential in the stream 
direction. The expression “linearized pressure” refers to the 
pressure obtained by neglecting all powers of the perturbation- 
velocity components above the first. 


VECTOR OPERATORS AND IDENTITIES 

Certain operators, which are closely associated with the 
linearized hyperbolic partial-differential equation of super- 
sonic flow (the two-dimensional wave equation), are added 
to the vector operators commonly used. The basic operators 
have been used previously in references 2 and 7. 

The gradient operator is defined by 


. d . d , d 

V = I r f“ J — b ^ x — 

dx dy dz 


The analogous hyperbolic gradient operator defined by 
Robinson in reference 2 may be expressed as 




dz 


The hyperbolic divergence of an arbitrary vector E is given by 

Vh-E 


Similarly, the hyperbolic curl of the vector E is given by 

VhXE 


The divergence of the gradient operator is sometimes denoted 

by 


V 2 


eV.V = 


+ ^+ 


d 2 

dz* 


The analogous divergence of the hyperbolic gradient operator 
is denotcdiby 


V 2 /r = V-VA= — p 2 



d 2 

dy 2 



The following identities are needed. Let E be a vector and 
\p and A be scalar functions of .r, ?/, and z. Then, 


Vi//-VA=VA-Vi/' 

(2a) 

V-fE=fV-E+E-V+ 

(2b) 

VX(VXj E)=V(V-E)-V 2 E 

(2c) 

V-(VXE )= 0 

(2d) 

V/i^-V \=VhA-Vt 

(2e) 

Vh.iE=iVhE+EVh^ 

(2f) 

VX (V/iX E)=Vh{V-E) — V 2 /i E 

(2g) 

V/iX (VXE)=V(Vh-E)-V 2 hE 

(2h) 

Vh-(VhXE) = 0 

(2i) 


These identities can be proved by direct expansion. 
The divergence theorem may be expressed as 


| E n da= j* V-Z? dv 


( 3 ) 
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where n is the normal unit vector to the element of area da. 
The vector n is expressed mathematically as 

n=iv l +jv 2 +kv 3 


teresting to note, however, that surfaces exist inside the 
volume of integration across which E n can be discontinuous 
while at the same time E n h remains continuous. It follows 
that for such a surface n and n h must satisfy the relation 


where p h v 2 , and v 3 are the direction cosines of the outward 
drawn normal to the element of area da. 

A theorem more general than the divergence theorem is 
given by (this theorem follows from the results of ref. 8 , p. 
87) 



viE x -\- Cov 2 Ey-\- C z v z E z )da 



d E z ~ d E v « 
dr + 2 dy + * 



dv 


where the subscripts x, y } and z refer to components of the 
vector E , and C Y , C 2 , and C z are arbitrary constants. Note 
that if C\ = C 2 =C 3 = 1 the preceding equation reduces to 
equation ( 3 ). If 

a = -/3 2 
C 2 =C z =l 


the preceding equation reduces to 


<j) (~ P ** E *+ rtE,)da=J (- /3 2 dv 


or 


where 


^)E-n h da — j Vh-E 


dv 


( 4 ) 


n h = — ifPvi +jv 2 +kv z 


If the divergence theorem as expressed by equation (3) 
is applied to a volume throughout which 

V*J?=0 

then the surface integral over the bounding surface is 

(jpE-n da = 0 


provided that no surfaces exist inside the volume of integra- 
tion across which the normal component of E is discontinuous. 
Similarly, if equation (4) is applied to a volume throughout 
which 

Vh-E= 0 


then the surface integral over the bounding surface is 


n n h = 0 


(5) 


Let Q( x ,y, z ) = ® be the equation of such a surface. 

n= y W+W+^? VQ 

and 




^ Qi 4 * Qy 1 + Qz 


jV hQ 


Then, 


where the subscripts indicate differentiation. Substituting 
the preceding expressions for n and n h into equation (5) yields 





( 6 ) 


Any solution of equation ( 6 ) set equal to zero is the equation 
of a surface across which X E may be discontinuous while 
X, r E remains continuous. The fact that the Mach cone 
from any arbitrary point satisfies equation ( 6 ) can be easily 
verified. The equation of the envelope of the Mach cones from 
an arbitrary line also satisfies equation ( 6 ) (ref. 9, p. 106). 

FINITE PART OF INTEGRALS WHICH ARISE 
IN STEADY SUPERSONIC FLOW 

In the following sections use is made of the concept of the 
finite part of an infinite integral. This concept was intro- 
duced by Hadamard (ref. 4) and has been used by a number 
of other investigators. The concept of the finite part is, 
however, sometimes confusing. This section was therefore 
included in an attempt to give a realistic picture of the 
finite-part concept and also to present the first steps of the 
derivation of the scalar and vector potentials. 

The concept of the finite part of double integrals as 
defined by Hadamard and used in this report is different 
from the concept of the finite part of double integrals as 
defined in reference 10 . The essential difference between 
these two definitions lies in the manner in which the singular 
points along the Mach cone are treated. 

In reference 3, page 183, a vector function is used in the 
derivation of the common scalar and vector potentials. 
The analogous vector function based on equation ( 1 ) is 

W=J { V 4 .- 0 V J t 

where 

R=J(x-tf-p(y- v y-p{z-i;y 


$ 


E n h da = 0 


The hyperbolic divergence of vector W with respect to 
variables £, 77 , and f is given by 


provided that there are no surfaces inside the volume of 
integration across which E-n h is discontinuous. It is in- 


V/*.JF=iv 2 /^=|(-0 2 


ay , ay , ay\ 
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The preceding equation indicates that the hyperbolic diver- 
gence of the vector W set equal to zero yields the partial- 
differential equation of linearized supersonic flow. A 
mathematical derivation of W can be obtained; however, 
for the purposes of this report such a derivation is not 
needed. 

The result of applying equation (4) to the vector W is 
V<t>-4X da= fl V 2 A</> dv (7) 


where T represents the area of integration when R=R'. 
The integral over the area T may be reduced to 

Jr W ( V</> ' n " + ^t) da ( 1 0) 

where r is given by 

rT 2 

Since R f is a constant, equation (10) can be written as 


When 0 satisfies equation (1) throughout the volume of 
integration, the right-hand side of equation (7) is zero; 
thus, 

(f>( ], V</> - da = 0 (8) 

wh en 


V 2 A0= 0 


Equation (7) is applied to a volume (denoted by v 0 ) 
enclosed in the forward Mach cone from the point (x,y, z). 
This volume is bounded by the surface given by R=R\ 
where R' is a small constant, and an arbitrary surface Si 
enclosed in the forward Mach cone from the point (x,y,z). 
A cross section of the region of integration is shown in 
figure 1. Note that this region is analogous to the region 
that is sometimes used in calculating the potential function 
satisfying Laplace’s equation (ref. 3, pp. 151-153). For 
regions such as the one shown in figure 1, equation (7) may 
be written as 

J T (jf Ti)' n * da+ f„G da 

= J ) { v>h« dv (9) 


Z 



S ? /r( V *' n * + V) 


da 


(id 


Equation (9) can now be written as 
P 2 <t> 


-m^r) ' u + J.,(s b) n ‘ *• 

= J 0 7 ? v2h4> dv ( 12 > 

If 0 is required to satisfy the linearized partial-differential 
equation of steady supersonic flow, then 


V 2 /t0 = () 


and equation (12) reduces to 

^ (V0.n,,+^) da+ | Q r V0-^»V ^)-n A da = 0 (13) 

If R r is made smaller and smaller the integrand of the integral 
over the area T in equation (13) remains finite except on the 
small area close to the point (, x,y,z ). In anticipation of tak- 
ing the limit of equation (13) as R' approaches zero, the small 
area close to the point (x,y,z) is removed from the area T. 
The area T is divided into two parts. One part is the area of 
T which is downstream of the surface given by 

t = X—6 

where e is small but larger than R f . This area is denoted by r. 
The remaining part of T (denoted by T') is the area of T 
which is upstream of the surface 


A cross section of the region of integration with T divided 
into r and T ' is shown in figure 2. Equation (13) can now be 
expressed as 


R' J( V *' n " + 



7?)'"" da+ 

j , G V ^ V ]?)"^ a =0 


where R' is smaller than e. 

Since 0 is continuous and therefore its values over r are 
approximately constant for small values of e, the integral 
over the area r can be written as 



U0-/1/, da- 


p 2 <t> 

R' 


l 


da 

r 


Figure 1. — Cress section of the region of integration used in connection 
with equation (9). 


(15) 
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Figure 2. — Cross section of the region of integration used in connection 
with equation (14). 

When the second integral of expression (15) is integrated, 
equation (14) becomes 

g? j V<t>-n h da + 2v — 2w<t>(x,y,z)+ 

g, da+ j* (t ^) n„ da= 0 

(16) 

If IV is made to approach zero, equation (16) applies even 
to the limit where R' is zero. 

The limit of equation (16) as IV approaches zero may be 
written as 

jnn Q*, J^V0-n„ da+~ Wi<t> ^,’ 1/ '~ ) —2Tr<t>(x,y,z) + 

R'fr'r </a + Js 1 (/? g) n " </a ] =0 (17) 

The integrands of t lie integrals over the areas r and T' are 
always finite and it can be shown that their first derivatives 
with respect to IV approach zero as IV approaches zero; there- 
fore, the product of 1 /IV and these integrals either approaches 
zero in at least the order of IV or approaches infinity as IV 
approaches zero. Thus it follows that the integrals over the 
areas r and T' have no finite terms remaining after the limit 
(/?'—» 0) has been taken. The sum of the terms of equation 
(17) must be zero; thus the singularities resulting from the 
integrals over the areas r and T' must cancel the singularities 
which arise from the integral over the area Si. 

From the preceding considerations it follows that one 
method of evaluating the finite part of infinite integrals of the 
type appearing in equation (17) is to evaluate the integral 
when IV is small but not zero and neglect t he terms multi- 


plied by powers of 1 /IV. Other infinite integrals sometimes 
arise, however, for which the finite part cannot be obtained 
by neglecting powers of 1 /IV. For example, if equation (17) 
is differentiated with respect to one of the variables ( x , y , or z) 
an equation containing the velocity component is obtained. 
In some cases, when the point (x.y,z) lies on the surface Sj 
the infinite terms are of the order (In IV) /IV and of the orders 
(l/IV) n . In these cases, the finite part of the infinite integrals 
can be obtained by evaluating the integrals when IV is small 
and neglecting the terms multiplied by powers of \/IV and 
(In R')/R'. 

The process of removing the infinite parts of an integral, 
however, has been derived by Hadamard (ref. 4, book III, 
ch. I). Hadamard used his methods of evaluating the finite 
part of integrals in finding solutions to certain hyperbolic 
equations including the linearized equation of steady super- 
sonic flow. Perhaps a fact worth noting is that t he integrals 
of equation (17) are double integrals and when the methods 
given by Hadamard are used the methods given for multiple 
integrals should be used. In the past, the singular points 
(points on the Mach cone where the derivative of 
(x—£) 2 —p 2 (y-v)) 2 —p 2 (z—{) 2 with respect to the variable 
of integration is zero) have caused some confusion; as 
Hadamard points out (ref. 4, p. 147), these singular points 
must be removed from the area of integration before the 
finite part is taken. Particular attention should be given to 
paragraph 92 of reference 4 since the special type of integrals 
discussed therein sometimes arises in dealing with planar 
problems. 

Robinson (ref. 2) has shown that when using Hadamard’s 
methods the order of integration may be changed without 
affecting the finite part and that it is permissible to differen- 
tiate under the integral sign of a multiple integral without 
considering the variable limits which lie along the boundary 
where the integrand is singular, provided that only the finite 
part is taken. Both Hadamard and Robinson have shown 
that in differentiating an improper integral which has an 
integrand that has a one-half power singularity along 
variable limits the variable limits may be neglected provided 
the finite part of the resulting integral is taken. 

The term “finite part” is somewhat misleading since the 
finite part of an integral can be infinite. In certain cases the 
integral is infinite even after the terms which approach 
infinity as R' approaches zero have been neglected. 

SCALAR POTENTIAL 

The preceding arguments show that the finite parts of 
equation (17) can be equated to zero; thus, 

— 2nt>(x,y,z)+ | ^ n h da = 0 

where the symbol / before the integral denotes that only the 
finite part is to be taken. The preceding equation may be 
solved for the value of the potential at the point (x,y,z); 
the result of this operation is given by 

<f>(x,y,z)=~ f ^ll) ' n,,<ia ( 18 ^ 

It should be remembered that surfaces can exist inside the 
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forward Mach cone from the point (x,y,z) across which 


1 1 
-^V0-0V-^ can 

• n„ 


be discontinuous 
remains continuous. 


and 


across which 


Equation (18) is an expression for the scalar potential at 
the point (x,y,z) in terms of the potential and its derivatives 
with respect to n n on the surface Si. A more general 
expression for the scalar potential than that given by 
equation (18) can be obtained. If within the volume 
enclosed by the forward Mach cone from the point ( x,y,z ) 
and the surface Si equation (1) is not satisfied and V 2 h 
operating upon 0 yields 


where S' 0 denotes the surface of discontinuity, and A 0 is 
the potential difference across the surface S 0 . The notation 
d/dft,/ is used to denote the operator 


1TZ +V/ 


b 


For the cases where no surface of discontinuity exists 
inside the volume V\ and 0 and V0 are zero on the surface 
S\ — So, equation (22) reduces to 


<l>(x,y,z)= 


JL f 

2 7T R 


dv 


From equation (19) 


(23) 


V 2 A0(f,iy,f) = /■({,? 7,f) 


(19) 


V 2 h<t>{x,y f z)=J(x,y,z) (24) 


then equation (10) becomes 

F X( v * ‘ n »+ 0 ~) da + 

= f f &gQdv (20) 

Provided that /(£,)/, f) is always finite, the right-hand side of 
equation (20) is finite; furthermore, the right-hand side of 
equation (20) remains finite as R' approaches zero. If in 
equation (20) R' is made to approach 'zero and onlvjthe finite 
parts of the integrals are retained, then the resulting expres- 
sion is 


Note that equation (23) is a solution of the partial-differential 
equation (24). 

For most problems in linearized supersonic flow, /(£,i 7 ,f) 
is zero and 0 is zero upstream of the disturbing body. For 
such problems, the surface Si—S Q can be taken to be located 
upstream of the disturbing body where 0 and V0 are zero. 
In this case, equation (22) reduces to 

J, (b * W S7 b) da (25) 

If the surface S 0 is confined to the f=0 plane, equation (25) 
reduces to equation (10) of reference 6. In this reference 
the boundary conditions for airfoils are discussed. 


4>(x 




dv 


( 21 ) 


where v x represents the volume when R ' is equal to zero. 
Equation (21) is equation (58) of reference 4 where fi 2 has 
been set equal to one. Note that the volume integral in 
equation (21) has the appearance of the integral for the 
potential resulting from a volume distribution of sources 
in an incompressible flow. 

The assumption has been made that 0 is continuous 
throughout the volume V\. It is also assumed that no 
surfaces exist inside i\ across which d<f>/dn h is discontinuous. 
If equation (21) is applied to a volume v u which has surfaces 
across which 0 and/or the derivative of 0 in the direction 
of n h is discontinuous, these surfaces of discontinuity can 
be removed from the volume of integration by allowing the 
arbitrary surface S\ to envelop them (see fig. 3). For vol- 
umes of integration where the surfaces of discontinuity 
have been removed in this manner, the scalar potential 
can be written as 


/ X 1 f 

J7T Jr t 




dr + — 


m 


i) i ‘+sL(s' # - 


i A d<t> - 
R dn h ' 

0V J , ) ■ n„ da 


COMPONENTS OF VECTOR FIELD 

Let F be a vector which is finite and integrable in a given 
volume (denoted by X) and is zero outside the volume X. 
To each point in the volume associate the vector 

G(x,y,z)= ( ^ dv (26) 

Jo 2 

where v 2 denotes the part of the volume X enclosed in the 
forward Mach cone from the point ( x } y,z ). 

From equations (24) and (26), it follows that each com- 
ponent of G satisfies the relation 

V 2 hG i (x,y,z) = —2rF i (x,y,z) (27) 

where the subscript i refers to any component of the vector 
field. 

Let 0o (:x,y,z) be a scalar and A(x,y,z) be a vector defined 
by the equations 

2*-^ =Vh-G= 'j F(H, n , n • VA ~ d v (28) 

and 

2^=VXG= / J | '"Ftf,v,r)XV-gd» (29) 

Equation (21i) indicates that 


( 22 ) 


VAX (VX G) =V(VA- G) — V 2 A G 


(30) 
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Figure 3. — A cross section of the region in the forward Mach cone 
from the point ( x,y,z ) showing the surface Si enveloping a surface 
of discontinuity So. 

Substituting the expressions for VXG, Vh»G, and V 2 hG 
as given by equations (27), (28), and (29), respectively, into 
equation (30) and solving for F yields 

F{x , y, z) = — Vxfy 0 (; x,y,z)+VhXA(x , y , z) (31) 

Since F is an arbitrary vector, equation (31) indicates that 
any finite integrable vector field can be expressed in terms 
of the gradient of a scalar and the hyperbolic curl of a vector. 
Equation (31) has the appearance of the Helmholtz theorem 
(ref. 3, p. 187); however, since \f/ Q and A are found by inte- 
gration only in the forward Mach cone from the point 
(x,y,z), equation (31) hardly seems to be a statement of 
the Helmholtz theorem as is commonly given. The result 
given by equation (31) was obtained by Robinson in 
reference 2. 

HYPERBOLIC VECTOR POTENTIAL 

Equation (31) indicates that the perturbation velocity 
vector can be divided into two parts. One part is the 
gradient of a scalar function, and the other is the hyperbolic 
curl of a vector function. The vector function is analogous 
to the common vector potential (ref. 3, pp. 104 and 188); 
therefore, the vector function is referred to as the hyperbolic 
vector potential. Thus, if q' is the total perturbation 
velocity vector, then 

q'=V(l)+VhXA (32) 

where 0 is the scalar potential and A is the hyperbolic* vector 
potential. The part of the velocity vector which is made up 
of the hyperbolic curl of the vector potential is denoted by q. 
By direct expansion it can be shown that 


Equation (33) indicates that the hyperbolic divergence of 
the perturbation velocity vector is zero. 

The vorticity vector is given by 

H=VXq' (34) 

Therefore, from equation (32), 

H=VX(VhXA) 

or 

H=Vh(V-A)-V 2 hA 

From equations (2d) and (29), the divergence of the hyper- 
bolic vector potential is zero; thus, 

H=—V 2 hA (35) 

Each component of equation (35) is a partial-differential 
equation of the form of equation (24); thus, from equation 
(23) each component of equation (35) has a solution given by 

A,(x,y,z)=-~ dv ( 36 ) 


where the subscript i refers to any component of the vector 
H. Since each component of A is given by equation (36), 


then 



g(€,n,» 

l R 


dv 


(37) 


The velocity vector resulting from the hyperbolic vector 
potential is therefore given by 


or 


, = v*x^-^j ri v»x f 

l ~2ir R 3 

~2 xjr, R 3 

V (x-QH,-(y-y)E 
2w J H " R 3 


dv 

(38) 

^dv 

(39a) 

Iz dv 

(39b) 

[z dv 

(39c) 


where the subscripts refer to the components of the vector 
H. The results given by equations (39) were obtained by 
Robinson in reference 2. 


VORTEX SHEETS 

If the vorticity is confined to a surface S 2 , equation (37) 
becomes 

da (40) 

Equation (40) is an expression for the hyperbolic vector 
potential resulting from a surface of vorticity. Note that 
if the vorticity is zero except on the surface S 2 , then equation 
(35) reduces to 

V 2 /*4 = 0 


Vh-q' = Vh-V<t>+Vh'(VhXA)=V 2 h<t>+Vh-(VhXA) = 0 (33) 


By removing the surface S 2 from the volume of integration 
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each component of A can be expressed as (from eq. (25)) 


. , , i 'f n . zAi . , b_i\ 

Mx,y,z)~2i J S 2 Vfl A £>V ‘d/o/ it) 


da ( 4 1 ) 


where the subscript i refers to any component of the vector 
A. Since each component of A is given by equation (41), 
then 

sih (42) 

Note that if A.4 is zero, equation (42) reduces to 

. \ W 1 , t>A . 

A(x,y,z)- 2v J ^ R A da 


(43) 


By comparing equations (40) and (43) it follows that, on the 
surface S 2 , 

//=a|4 (44) 

£>ri n 

Equation (44) indicates that across a surface of vorticity 
the derivative of the hyperbolic vector potential in the di- 
rection of rtn is discontinuous. Thus, a lifting surface can 
be represented by a continuous hyperbolic vector potential, 
while it can be shown that a thickness effect can be repre- 
sented by a discontinuous hyperbolic vector potential. 
Note the contrast with the scalar potential, which uses a 
continuous potential to represent a thickness effect and a 
discontinuous potential to represent a lifting surface. 

FURTHER DEVELOPMENT OF SCALAR POTENTIAL 

The scalar potential can be expressed in forms other than 
th ose already presented. Equation ( 8 ) is applied to the 
region bounded by the arbitrary surface S\, the forward 
Mach cone from the point (x,y,z), and a second arbitrary 
surface £3 enclosed in the forward Mach cone from the point 
(x,y,z) and upstream of the surface S,. A cross section of 
such a region is shown in figure 4. The result of applying 
equation ( 8 ) to this region is 


Ts.G 


d</>' 


-</>' 


r 


R bn* v dc/ It ) ' la + J s,( 


1 d<f>' 


-<j) 


, _J> J\ 

RJ 


dZ* i>n»* Rr a ° 


Equations (18) and (46) can be combined to yield 




da =0 

(IS) 


provided that </>' is a solution of equation ( 1 ). Note that the 
scalar potential as given by equation (18) is independent of 
</>' so that </>' is arbitrary so long as it satisfies equation ( 1 ) 
throughout the proper volume. 

If for a finite distance upstream 0 ' is zero and remains zero 
for greater distances upstream, the surface may be chosen 
in this region so that the integral over S s in equation (45) is 
zero; thus, 

J_ r / 1 d 1 \ 

27 r J sXR 


(46) 



Figure 4. — Cross section of the region of integration used in connection 
with the potential function <f>' . 

The only restrictions placed on <£' at this point are that it 
satisfy equation ( 1 ) and be zero at a finite distance upstream. 

def) C )<// 


In many cases <j > ' mav be chosen so that r- — ■- ---* or <f> — 

c )n h dn h * 

is zero; therefore, in these cases, cf> can be expressed as 


or 


\ 1 r 1 / c )<t> c)(f>' \ , 

<Hx,y,z)~ j sR da 

*(x,y,z)=-± 'J s &-*')*- 


d 1 
R 


(48a) 


(48b) 


(47) 


Equations (48) are quite useful; however, remember that 
they apply only when <f>' can be chosen so that 0 ' does not 
violate any of its restrictions. 

Note that equations (48) can be applied to problems where 
either 0 or 50/dn,, is given on the surface S 1 . The application 
of these equations to most nonplanar problems of either type, 
however, lead to quite unwieldy integral equations. 

APPLICATIONS 

PLANAR PROBLEMS 

Many problems in linearized supersonic flow deal with the 
surface of discontinuity confined to a plane surface parallel 
to the j:-axis. In this section a general discussion of this 
type of problem is given. The coordinates are located so 
that the surface of discontinuity is in the plane. 

The scalar potential at an arbitrary point (x,y,z) above the 
plane is (from eq. (47)) 

j^( 0 2 — 0 /) j^+( 0 — 0 ') d£ drj 


( 49 ) 


In this case, the surface S 1 is the f = 0 plane. 
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If <{>' (x,?y,0) is chosen equal to <f)(x,y y O) the potential becomes 
symmetric with respect to the f=0 plane. Thus, for z=a 


and 


4>(.x,y } a) = 4>' (x,y,—a) 
<t>z(x,y,a)=—<t>/(x,y,—a) 


For this case, equation (49) reduces to 


dt dn (50) 

Equation (50) was given by Puckett in reference 5. 

If <t>'(x,y, 0) is chosen equal to — <t>{x,y, 0) the potential 
becomes antisymmetric with respect to the f=0 plane. In 
this case, 

<#>(•»•', y, a) = — <*>' (x,y, -a) 

and 

<t> z (x,y,a)=<l>/(x,y,—a) 

Thus, equation (49) is reduced to 


0Or,?/,2)= 


7 T J S, R 3 


dt; dn) 


(51) 


Note that for surface S { not confined to a plane parallel to the 
$-axis, a choice of <j>'(x,y,z) at the surface S\ to equal <t>(x y y,z ) 

at Si does not cause at the surface S x to equal 

dn h 1 

_ d0 (x,j/, 2 ) ^ Similarly, choosing <t>'(x y y,z) at the surface 

on h 

Si equal to —ct>(x,y,z) at S x does not cause ^ — at S\ to 

OW/t 

equal at 8 t . 

Provided the discontinuities are restricted to the f=0 
plane, the scalar potential can also be expressed as follows 
(from eq. (18)): 


</>(«, 


2 TV 



$*(£> V >0) 

R 





(52) 


for positive z. A comparison of equations (50) and (51) with 
equation (52) shows that the two terms of the integrand of 
equation (52) contribute equal amounts to the potential at 
any point ( x,y,z ). 

Since the terms of the integrand of equation (52) contribute 
equal amounts to the potential at the point (x,y,z) as 0 
approaches zero, equation (52) must reduce to 


<Kx,y, 0 )= 


_L f 4>z(U y ,0) j ,(t>(x,y,0) 
2 ttJs 1 R <l - (71 + 2 


The preceding equation can also be obtained by examining 
the limit of equation (52) as 2 approaches zero. If this 
procedure is done the entire contribution of the second term 
of the integrand of equation (52) is found to come from the 
point at the apex of the hyperbola formed by the intersection 
of the Mach cone from the point (x,y,z) and the f=0 plane. 
Note that if the integration is performed first with respect to 


7] then, when the methods of Hadamard are used, the point at 
the apex of the hyperbola is a singular point and must be 
removed from the area of integration by a process such as is 
given in reference 4, page 147. 

If 0 2 is prescribed over the f=0 plane, then the potential 
is given uniquely by equation (50). Similarly, if the potential 
is prescribed over the f=0 plane, then the derivative of 0 with 
respect to 2 is determined over the 2=0 plane. This result 
follows from equation (51) since prescribing 0 over the 2=0 
plane determines the potentials in the space above the 2=0 
plane; therefore, it also determines 0, in the space above the 
2 = 0 plane and the limit of 0, as 2 approaches zero from the 
positive direction. 

The question that arises is whether 4>(x } y,z) is uniquely 
determined in the space above the 2=0 plane if </> is prescribed 
over certain areas of the 2=0 plane and <t> 2 is prescribed over 
the remaining areas. If the assumption is made that <f> is 
not determined uniquely, then at least two potential functions 
satisfy the condition that either 0 or 0 2 is prescribed in all 
regions on the 2=0 plane and that 0 is identically zero 
upstream of a given point. Let 4>, and # 2 denote two poten- 
tial functions which satisfy the same boundary conditions, 
and let <t> 0 denote the potential function formed by taking the 
difference between </> j and <f> 2 . Mathematically, the potential 
function <£ 0 is given by 


<t> 0 (x,y,z)=<l> l (x,y,z) — <t> 2 tx,y,z) (53) 

Since <t> { and <p 2 have the same values in certain regions in the 
2=0 plane then <l> 0 is zero in these regions. Similarly, since 
dfa/dz and d<£ 2 /d 2 have the same values in the remaining 
regions of the 2=0 plane, then d</> 0 /d 2 is zero in these remain- 
ing regions. The potential function <j> 0 has the boundary 
conditions that either 4> 0 or d<t> 0 /dz is zero in all regions of the 
2=0 plane and that 4> n is identically zero upstream of a given 
point. 

Consider the case where all the boundaries between the 
regions are supersonic. (The slopes of the boundaries are 
such that the component of the free stream perpendicular 
to the boundary is always supersonic.) The potential 
function 0 O can be evaluated by use of equation (50) or (51) 
for points in areas which are far enough upstream to be 
affected only by a region where 0 or 0 : is prescribed. For all 
points in these areas 0 O is zero as indicated by equation (50) 
or (51). It follows from equations (50) and (51) that 0 O is 
also zero inside the volume above the 2=0 plane, which is 
affected by these areas alone. Thus, the volume where 0 O 
is identically zero has been moved downstream. The same 
argument can be repeated until the complete 2=0 plane has 
been covered. The preceding arguments cannot be applied 
to cases where the regions have subsonic boundaries; however, 
if it is permissible to distort the boundaries within a strip of 
infinitesimal width these subsonic edges can be converted 
into supersonic edges by replacing every element of the 
subsonic boundaries by a broken line made up of supersonic 
segments. Such a procedure is illustrated in figure 5. If 
the assumption is made that the subsonic boundaries may be 
distorted an infinitesimal amount, then 0 O is zero over the 
2=0 plane and also in the space above the 2 = 0 plane. 
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(a) Subsonic element within (b) Broken line replacing 
infinitesimal strip. line element. 

Figure 5 . — An illustration of a method of replacing a subsonic line 
element with a broken line made up of supersonic segments. 

Equation (53) now reduces to 

<t> i(x,y,z)=<t>i(x,y,z) 

Since <t>\ and 0 2 are any two potential functions with the same 
values in certain regions of the z= 0 plane and with the same 
partial derivatives with respect to 3 in the remaining regions, 
proof has been given that only one potential function exists 
for which the potential is prescribed over certain areas in the 
z= 0 plane and the partial derivative with respect to z is 
prescribed over the remaining areas. 

The boundary conditions for a zero- thickness lifting airfoil 
with a given local angle-of-attack distribution are not of the 
type discussed in the preceding paragraph. The conditions 
prescribed in the z=0 plane for this type of problem are: 
The potential function 0 is identically zero upstream of the 
airfoil; <t>(x,y y 0) is zero except on the plan form or in the 
wake; the partial derivative of the potential with respect to 
z , 0,, is given on the plan form; and <t> x (x } y,0) is zero in the 
wake. The preceding boundary conditions do not specify 
that 0 or 07 be prescribed in all regions on the z=0 plane 
since not 0 but <t> x is given in the wake. For airfoils which 
have trailing edges which are always supersonic, the require- 
ment that 0 be continuous in the stream direction necessi- 
tates the potential in the wake to have the value of the 
potential at the trailing edge of the airfoil. In this case, the 
potential function is uniquely determined. For airfoils 
which have subsonic trailing edges the Kutta-Joukowski 
condition is generally applied to the trailing edges to deter- 
mine 0 uniquely. If the assumption is made that the trail- 
ing edge can be distorted within a strip of infinitesimal width, 
then the requirement that 0 be continuous in the stream 
direction can be used to determine 0 uniquely. If the as- 
sumption is made that the subsonic trailing edge is distorted 
within the infinitesimal strip so that each segment of each 
line element of the trailing edge is always supersonic (see 
fig. 6), then 0 is determined uniquely. It is well-known 
that for airfoils with subsonic trailing edges there are an 
infinite number of solutions which satisfy the boundary 
conditions as stated at the beginning of this paragraph. The 
preceding arguments however prove that there is but one 
solution for an airfoil which has had its subsonic edge re- 
placed by broken lines which are always supersonic. Note 
that it has not been proved that the solution obtained by 
distorting the subsonic trailing edges corresponds to the 




(a) Wing with subsonic edge. (b) Line element at a break 

in a subsonic trailing edge. 

(c) Possible distortion of a subsonic trailing edge at a break in 

that edge. 

Figure fi. — Methods of distorting a subsonic trailing edge to determine 
the potential function uniquely. 

solution satisfying the Ivutta-Joukowski condition, nor has 
it been proved that the solution of the distorted trailing edge 
is independent of the manner of distortion. 

NONPLANAR PROBLEMS 

The scalar potential resulting from the disturbances caused 
by a nonplanar body can be found from equation (18) 
provided that both 0 and 50/dft/, are known on some surface 
S\. Unfortunately, 0 and d<t>/dn h are not generally known 
on a surface which fills the requirement of the surface Si; 
therefore, equation (18) appears to have little value in the 
calculation of the potential functions for nonplanar systems 
in general. Certain properties of equation (18) are, however, 
worth investigating. 

The problem of evaluating the potential on the upper sur- 
faces of a long rectangular body is discussed. The assump- 
tion is made that the body extends upstream to infinity and 
that the sides are parallel to the free-stream direction except 
for small local variations which cause small disturbances in 
the stream. Figure 7 (a) shows the forward Mach cone from 
a point on the upper surface of such a rectangular body. 
This figure also shows that there is a certain part of the sur- 
face of the rectangular body in the forward Mach cone from 
the point ( x,y,z ) that cannot possibly affect the potential at 
the point (x,y,z). If the surface S x in equation (18) is 
taken to be the surface of the rectangular body, then equation 
(18) indicates that the values of 0 and 50/5/i^ in the region 
which cannot possibly affect the potential at the point 
(x,y,z) should be used in evaluating the potential at the 
point (x,y,z). The only possible explanation of this con- 
sideration is that the integral of 0 and d(t>/dn n caused by the 
disturbances in the “blind spot” add to zero. This con- 
sideration can be shown mathematically as follows. Let 0 O 
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denote the potential function resulting from the disturbances 
inside the blind spot. From equation (46), it follows that 


T /1_ djAo 
J s\R bn,* 




da = 0 


(54) 


Equation (54) indicates that the potential at the point 
(x,y,z) can be evaluated by applying equation (18) to the 
surface of the rectangular body regardless of the blind spots. 
The same argument holds for other bodies with blind spots. 

The preceding arguments can be clarified by a simple 
illustration of the effect of a blind spot. Consider an infinite 
rectangular body such as shown in figure 7 where the only 
disturbances are caused by a small deflected area with a 
constant slope a with respect to the free-stream direction 
located on the lower surface of the body. The leading edge 
of the deflected area is chosen perpendicular to the free- 
stream direction so that the potential in the region not 
affected by the vertical sides is of a two-dimensional nature. 
Figure 7 (b) illustrates such a disturbing surface. 

The disturbance potential in the two-dimensional region is 
given by 

, , „\_(rV[x—x 2 +0(z — zJ] 



(c) Side view of rectangular body and the forward Mach cone from the 

point (x,y,z). 

Figure 7. — Concluded. 


A point on the upper surface which has only the two- 
dimensional flow in its forward Mach cone is illustrated in 
figure 7 (c). The disturbance potential for this point is 
(from eq. (18)) given by 


</>(x, 2 /,z)=— 

Z 7T 



P 2 (tV(£—x 2 ) (z — 

R? 



where the lower surface of the body lies in the z=z 2 plane 
and the leading edge of the deflected area is in the x=x 2 
plane. 



z 



(a) The forward Mach cone from a point on the upper surface of a 
rectangular body. 

Figure 7. — Rectangular body parallel to free-stream direction. 


Deflected area 




5 



(b) 


(b) Bottom view of rectangular body with deflected area on the lower 

surface. 

Figure 7. — Continued. 


Upon performing the indicated integrations the preceding 
expression becomes 

j <rV[x—X 2 —P(z — Z 2 )] , aV[x — X 2 — 18(2 — Z 2 )] 

4>{x,y,*)= 2 ^ 20 - 

which reduces to 

<t>{x,y,z) = 0 

This result is a demonstration that the disturbances in blind 
spots do not contribute to the potential. 

The scalar potential resulting from the disturbance pro- 
duced by a nonplanar body can also be obtained by use of 
equations (47) and (48) provided that the necessary values 
of </>, b(f)/bn h , <£', and b4>'/bn h are known. In dealing with 
planar bodies <£' could be chosen so that 

cf ) — <t>' = 0 

and thus equation (47) is reduced to equation (50). Simi- 
larly, </>' could also be chosen so that 

b<j> b<t>' 
dn h ~dn h * 

and thus equation (47) is reduced to equation (51). Un- 
fortunately, for nonplanar bodies, choosing </>' equal to </> 
does not make b<t>'/bn f * known as was the case for planar 
problems and, similarly, choosing b())'/bn h * equal to — b<t>/bn h 
does not make </>' known. Certain problems exist in which 
</>' can be chosen so that 0 can be written as a simple integral. 

INTERSECTING PLANES 

Many problems concerning nonplanar bodies deal with 
disturbances produced by two intersecting planes parallel 
to the free-stream direction. In this section, methods of 
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solutions for two pianos intersecting at various angles are 
given. The component of velocity normal to the surface is 
assumed to be known. 

Perhaps the simplest case of two intersecting planes occurs 
when the planes intersect at right angles. It is desired to 
find the potential in space resulting from the disturbances 
produced by the two intersecting planes. This type of 
problem could represent an isolated cruciform tail with 
supersonic leading edges undergoing various motions. 
Problems of this type have been solved in references 1 1 and 12. 
The axes are chosen so that y = 0 and 2 = 0 are the disturbing 
planes (see fig. 8). When y and z are positive, equation 
(18) becomes 


tf>(z,y,z)= — — 
27 r 


f 1 fy>(s,o + ,r) 

J «3 R &V 




R'-' 


da — 


J_ C ld0(£,„.O + ) 
2 ir J s 4 R d { 


'-g'LV 


da 


(55) 


The surface S\ has been taken to be the disturbing surface; 
thus, S 3 is the part of the y= 0 plane (z positive) bounded 
by the 2=0 line in the y= 0 plane and the forward Mach 
cone from the point (x, y, 2 ). Similarly, S 4 is the part of the 
z=() plane (y positive) bounded by the y=0 line in the 2 = 0 
plane and the trace of the forward Mach cone from the point 
(2, y. 2 ) (see fig. 9). 

11 0 + ,f) , c></)($,^7, 0 + ) 

1 he assumption is made that ^ — and — are 

known and that </>(£, 0 + ,f) and </>(£, 77, 0 + ) are unknown. The 
integrals containing <M£,0 + ,f) and 0(£,77,O + ) can be eliminated by 
several applications of equation (46). Equation (46) is 
applied to the volume on the left-hand side of the y = 0 plane 



enclosed by the forward Mach cone from the point (x,y,z), 
the y=0 plane, the 2=0 plane, and an arbitrary surface 
upstream of the disturbance (see fig. 10). The result of 
applying equation (46) to this volume is 


1 [• 1 d</>u^o+) , isvf 

~2tJs 5 R dr 2 ttJs, 

a. 


<*>'($, 0 +,r) 

, R 3 


da-\- 




rfo +£vr <(L>±:>n da 

® + 2 t r.)s 3 R 3 


= 0 (56) 


The surface S\ has been taken to be the y— 0 plane (2 posi- 
tive) and the 2=0 plane (y negative); thus, S 5 is the part of 
the 2=0 plane in the forward Mach cone from the point 
(j\y,2) (see fig. 10). Adding equations (55) and (56) yields 


J(rv .) 1 r ip*({,o+,r) -,m. 

"2^J,3«L ^ V 

^>1 \ h [ ^’ 1 0 ~ 0 f )] ■ da ~ 


da — 


y, 0 + ) 

R 

pz'C i-o + ) . L f ar 

2 ir J s 4 R 3 0(1 2 xjs. R 




2 

r*'r 

2tt J .< 


da — 


4 >'U, 0 + ,f) 

. 


r/a 


(57 1 



Figure 8. — Two disturbing surfaces intersecting at right angles. 


Figure 9. — Regions of integration for equation (55). 


A VECTOR STUDY OF LINEARIZED SUPERSONIC FLOW APPLICATIONS TO NONPLANAR PROBLEMS 


13 


The potential function </>' (£,*;,£) is chosen so that 
</>'(£, — a, f) = <!>(£, a, {) 
where a is positive. In this case, 

a+'foo- r) = a»foo + ,r) 

dr/ dr/ 

and equation (57) reduces to 


I f 

1 o<t>a, o+.r) 

d a L 

f 1 0<t>(£,T), 0 + ) 

da- 

TrJ S3 

R Or] 

da 2ir, 

)s 4 R df 

0*2' 

r <ks,v, o + ) 

c la — — I 

r 1 &*'(€, ij, o + ) 

da- 


1 «4 R 3 

2 rrj 

sj{ df 

Fz j 

2ir J 

'• 0 + ) 
« 5 R 3 

da 


rn 


Since <£'({, v,t) is related to <!>(£, rj,{) the only unknown in the 
preceding equation is <#>(£, *7, 0 + )- The region of integration 
S 5 becomes the part of the 2 = 0 plane (2 negative) in the 
forward Mach cone from the point ( x,y,z ) obtained by re- 
flecting the disturbing surface in the 2=0 plane (z positive) 
through the 2=0 plane (see fig. 10). 

The problem being considered is one in which the normal 
derivatives of the potential function are known on two planes 
parallel to the #-axis and are intersecting at right angles. 
The point ( x,y,z ) has been restricted to positive values of 
y and 2 . For the present, consider the problem of finding 
the potential above the 2=0 plane when the derivative of 



the potential function with respect to 2 is known on the 2=0 
plane, the derivative of the potential function with respect 
to y is discontinuous across the y = 0 plane, and the potential 
is zero a finite distance upstream of a given point. From 
equation (22), the potential for this problem is given by 


1 f 1 A d </>(£, 0,f) , 

^ x ’y^=- 2 ^j So R A da ~ 


-f 1 

2t J st I) 


r dr 


da 


2 t r J sT 


R 3 


da (59) 


where S 0 is the part of the y= 0 plane above the 2=0 line 
inside the forward Mach cone from the point ( x,y,z ), and S { 
is the part of the 2=0 plane inside the forward Mach cone 
from the point ( x,y,z ). For positive values of y and 2 , 
equation (59) reduces to equation (58) if </>'({,?/, 0 + ) is assumed 
to be the true value of 0 (f, 7 /,O + ) when 77 is negative. Since 
in the original problem </>(£,??, r) was not defined when 77 was 
negative, nothing is violated if it is now defined as being 
</>'(£, 77 , r) ill tile region where 77 is negative and r is positive. 
The problem in which the normal derivative of the potential 
function is known on two planes parallel to the £-axis and 
intersecting at right angles has, therefore, been changed to 
the problem in which the derivative of the potential has a 
known discontinuity across the y = 0 plane (2 positive) and 
the normal derivative of the potential is known on the 2=0 
plane. Note that the potential function still remains unde- 
fined below the 2=0 plane. 

Since the potential function is undefined below the 2=0 
plane, it can now be defined so that the resulting potential 
function is symmetric with respect to the 2=0 plane. De- 
fining the potential below the 2=0 plane so that 


yields the desired symmetry. The result of applying equa- 
tion (46) to the region below the 2=0 plane inside the forward 
Mach cone from the point ( x,y,z ) is 


1_ 
2 7 T 

J_ 

2 7 T 


S t R Os 2ir 

* 1 d4>(£,7 ? ,0+) . Pz 

s S R Os 


Y <£(£> 0 + ) 

J «4 R 3 

T ^ Q + ) 
.U R* 


(Id “f - 
da-\- 


ir i a<Ks,o + ,f) 

T J S, R Or] 


da — 0 


(60) 


where S ( , is the part of the y—0 plane (2 negative) inside the 
forward Mach cone from the point (x,y,z). In applying 
equation (46) to the region below the 2=0 plane, the surface 
Si had to be folded over the part of the y = 0 plane (2 nega- 
tive) across which 04>/0r] was discontinuous in order to be 
able to apply equation (46) to this region. Equations (58) 
and (60) can be combined to yield 




( 61 ) 
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Equation (61) contains only integrals of known expressions, 
and it is, therefore, the solution to the problem of two planes 
intersecting at right angles parallel to the z-axis where the 
normal derivative of the potential function is known on both 
planes. Figure 1 1 shows the cross section of the distribution 
of velocity normal to the surfaces for a problem as repre- 
sented by equation ( 55 ) and its solution as given by equation 

(61). Note that (77 negative) in S 5 is the reflection 

of * in Si (77 positive) across the 77 = 0 plane and 

of 

that ^ s,() in S 6 (f negative) is the reflection of 

Orj Or] 

in S$ (f positive) across the f =0 plane. This condition 
suggests that the result given by equation (61) could also be 
obtained by utilizing the concept of reflecting surfaces. 

The mathematical derivation required for finding solutions 
to problems consisting of two planes parallel to the z-axis 
intersecting at various angles can be reduced by making use 
of the concept of the reflecting surfaces. For this reason, the 
result given by equation (61) is obtained by use of reflecting 


z 



z 



(a) Original problem (equation (55)). 

(b) Solution to problem as given by equation (61). 

Figure 11. — A cross section of the distribution of the velocity com- 
ponent normal to the 2 = 0 and the y = 0 planes represented by 
equations (55) and (61). 


surfaces. The potential function can be separated into two 
parts, 0 1 and fa, satisfying the following boundary conditions 
on the disturbing surfaces: 

d0i(£»q>O + ) 7 ,0*) 

df df 

d V 

d <&>(£, *7,0+) 

df 

dr] dr] 


A cross section of these boundary conditions is shown in 
figure 12. Only the potential function fa is treated in detail 
since the boundary conditions for fa and fa are of the same 
type. The normal derivative of fa is zero on the 77 = 0 
plane; thus, the rj = 0 plane can be considered as a reflecting 
plane. The potential function fa is, therefore, the potential 

resulting from a distribution of which is symmetric 

with respect to the 77 = 0 line and has the value of ) 


when rj is positive. Figure 13 illustrates such a distribution. 
The problem of finding fa has been reduced to a planar 
problem which can be solved by use of equation (50). 

Equation (50) was obtained by defining the potential 
below the 2=0 plane so that the total potential function was 
symmetric with respect to the 2=0 plane. This result 
caused the derivative of the potential function with respect 
to 2 to be antisymmetric with respect to the 2=0 plane. 
Figure 14 illustrates the distribution of the normal derivative 
of the function across the 2=0 plane. The problem of 
evaluating the potential function fa has been reduced to a 
planar problem. Similarly, the problem of evaluating the 
potential function fa can be reduced to a planar problem. 
Figure 15 illustrates such a procedure. The original poten- 
tial function is the sum of fa and fa. Equation (61) follows 
from the preceding results for 0, and fa. The addition of 
0! and fa is illustrated by figure 16. 

The concept of reflecting surfaces is now utilized to find the 
potential resulting from two disturbing surfaces parallel to 
the z-axis and intersecting at an angle of 45°. The axis is 
chosen so that the z-axis lies along the intersection of the 
disturbing surfaces and one of the disturbing surfaces lies 
in the 2=0 plane (see fig. 17). The potential function 0 is 
divided into two parts, fa and fa. The boundary conditions 
on 0i and fa are similar to the corresponding potential func- 
tions used for the disturbing surfaces intersecting at 90°. 
Figure 18 illustrates the boundary conditions for fa and fa. 
The surfaces on which the normal derivative of 0 1 is zero 
can be considered as a reflecting surface. This consideration 
leads to the same distribution of the normal derivative of 
0i on the 77 = 0 plane as is given on the f = 0 plane. Figure 
19 illustrates such a distribution. The problem of finding 
0! for two disturbing surfaces intersecting at 45° has been 
reduced to a problem of two surfaces intersecting at 90°. 
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Figure 12. — A cross section of the distribution of velocity normal to the disturbing surfaces for the potential functions <t> } 4> u and fc- 
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Figure 13. — The reduction of fa to a planar problem. 




Figure 14. — The normal derivative of the potential function across the 2=0 plane obtained by applying equation (50) to a planar problem. 


The solution of fa can be obtained from equation (61). 
Figure 20 shows the surfaces across which the normal deriva- 
tive of fa is discontinuous. Since fa and fa have the same 
type of boundary conditions, then fa lias a solution as illu- 
trated in figure 21. The original potential function <j> is the 
sum of fa and fa ; therefore, </> can be found by considering 
surfaces of discontinuity as illustrated in figure 22. The 
potential function </> can be evaluated by use of equation 
(22), because no surfaces across which </> is discontinuous 
exist and the values of Aid^/drifi) are known across all sur- 
faces of discontinuity. 

Another simple case of two disturbing surfaces parallel to 
the z-axis occurs when the surfaces intersect at an angle of 
60°. The potential function <£ is divided into two parts, 
fa and fa. The boundary conditions on fa and fa are similar 
to the corresponding functions used previously. Figure 23 
shows a cross section of these boundary conditions. By use 
of a reflecting surface, the function fa can be represented by 
the boundary conditions as shown in figure 24. The func- 
tion fa is undefined for 240° of the total angle around the 
£-axis. The function fa is defined as shown in figure 25. 


Since no surfaces exist across which fa is discontinuous, the 
function fa can be evaluated by using equation (22). Simi- 
larly, fa can be defined as shown in figure 26. The sum of 
fa and fa is illustrated by figure 27. The potential function 
(j> can be found by using equation (22). 

In the preceding paragraphs, methods have been found 
for determining the potential resulting from two plane dis- 
turbing surfaces parallel to the stream direction intersecting 
at certain angles. The same method can be used to find 
methods for determining the potential resulting from two 
plane disturbing surfaces intersecting at various other angles. 

ROLLING TAILS WITH MULTIPLE RECTANGULAR FINS 

The methods derived in the preceding section are used to 
find the surface velocity potential, the pressure distribution, 
and the damping in roll of rolling tails consisting of four, 
six, and eight rectangular fins. For comparison, these same 
quantities are also presented for the planar tail configura- 
tions consisting of one and two rectangular fins. An illus- 
tration of the tails treated is shown in figure 28. The 
analysis is limited to tail configurations having surfaces of 


276818—54 
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Figure 15. — An illustration of the reduction of <f> 2 to a planar problem. 






ggnmnM: 






Figure 16. — An illustration of the addition of <t>\ and <f> 2 . 


vanishingly small thickness and of zero camber. The 
investigation is also limited to the range of Mach numbers 
for which the region of interference between the adjacent 
fins does not affect the fin tips (see fig. 29). 

TAIL CONSISTING OF ONE FIN 

The pressure distribution and the velocity potential on 
the surface of rolling tails made up of one and two rectangular 
fins can be obtained from the results of reference 13. The 
pressure and potential for the tail consisting of only one fin 
can be found by transforming the axis of roll of the tail 
consisting of two rectangular fins. 

The tail consisting of one fin is divided into regions as 
shown in figure 29 (a). The velocity potential on the sur- 


face facing the negative ?/-direction is given for the 
regions as follows: 

For region I, 


2v 

7 r 



K¥~»)VG-)(-+H). 


For region II, 


</>(£, 2) 


pxz 


various 

(62a) 

(62b) 


A VECTOR STUDY OF LINEARIZED SUPERSONIC FLOW APPLICATIONS TO NONPLANAR PROBLEMS 


17 



igure 17. — Position of coordinate axes for disturbing surfaces 
intersecting at 45°. 


For region III, 

* (r ' z) ~ir [f sin "Vf (f + 0 V" Gr *)] (62c) 

For region IV (note that the potential in region IV is the 
potential in region I plus the potential in region 111 minus 
the potential in region II), 


4>( x > z)=— 

7T 


xz 


L 0 


sin 


i 3z . xz . 

Vt+T sm ' 



Kt+Ov-G-*)- 


Kt— < 0V(i-)(- 


b x\ irxz 

(62d) 


The pressure-difference coefficient is given for the various 
regions as follows: 





Figure 18. — A cross section of the velocity distribution on the disturbing surfaces for the functions <f>i and </> 2 for the 

disturbing surfaces intersecting at 45°. 



Figure 19. — Reflection of the normal derivative of <j > i on the y = 0 plane. 
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Figure 20. — A cross section of the surfaces of discontinuity which can be used to evaluate the potential function <f>i for 

two surfaces intersecting at 45°. 


For region I, 

A C = 

* ' tt(3V 

For region II, 


For region III, 


— 

l/b \ 



2 sill" 1 *! 

/ \ * 

L } 

)l * 


A C = *21 

' pv 




For region IV, 
AC =J&- 

" TpV 


Z Sill 


- 0 -*) 


+ 2si,l " , Vf + 


V*(H-V(F)H+0-¥ j 


(63a) 

(63b) 

(63c) 


(63 d) 


TAIL CONSISTING OF TWO FINS 

The tail consisting of two fins (fig. 29 (b)) has the same 
potential and pressure distribution as a rectangular rolling 
wing and can, therefore, be obtained from the results of 
reference 13. For each tail consisting of two fins divided 
into regions as shown in figure 29 (b) the pressure and poten- 
tial in regions I and II are the same as the pressure and po- 
tentials in the corresponding regions for tails consisting of 
one fin. 

TAIL CONSISTING OF FOUR FINS 

Each fin of the tail consisting of four fins is divided into 
regions as shown in figure 29 (c). The pressure and poten- 
tials in regions I and II are the same as the pressure and 
potentials in the corresponding regions for tails consisting of 
either one or two fins. The regions III and IV are affected 
by the interaction between adjacent fins. The potential in 
region IV is made up of a combination of the potentials of 
regions I, II, and III. Thus, the only real problem is the 
determination of the potential in region III. 

The potential in region III is not affected by the tip and 
is, therefore, the same potential as would be obtained if the 
fins were infinitely long. With the coordinate axes chosen 
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as shown in figure 28 (c) the point ( x,y,z ) is restricted to 
values of y which are negative while the values of z are 
restricted to positive values. Note that for a tail with 
finitely long fins, the potential at a point ( x,y,z ) in the 
region of interaction is independent of the disturbances 
produced at points located so that their projection on the 
yz - plane does not lie in the second quadrant. The general 
method previously derived for finding the potential resulting 
from two plane surfaces intersecting at right angles can thus 
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be used to find the potential in the part of the region of 
interaction which is not affected by the tip. 

The velocity component normal to the fin in the z=0 
plane is given by 

<t>z(x,y,0 + ) = -py 

and the velocity component normal to the fin in the y= 0 
plane is given by 

<t> y (xfi-,z)=pz 




Figure 21. — A cross section of the surfaces of discontinuity which can be used to evaluate the potential function fc for 

two planes intersecting at 45°. 




z 



Figure 22. — The addition of fa and fa to obtain the potential function $ for two surfaces intersecting at 45°. 
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Figure 23. — A cross section of the boundary conditions of the functions <f>\ and fa for two disturbing surfaces intersecting at 60°. 




Figure 24. — The reflection of <£i through one surface. 




Figure 25. — A method of defining fa so as to eliminate discontinuities in the potential function. 
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Figure 27. — The addition of <f>i and <f> 2 . 
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(a) One fin. (b) Two fins, 

(c) Four fins. (d) Six fins. 

(e) Eight fins. 

Figure 28. — Types of tails treated. 


Figure 30 illustrates this typo of normal-velocity distribution, 
and figure 31 illustrates a cross section of the surfaces of 
normal-velocity discontinuity, which previous results show 
can be used to obtain the potential in the part of the region 
of interaction which is not affected by the tips. 

Note that in figure 32 the discontinuity in the normal 
velocity across the 2=0 plane is the same type as the dis- 
continuity in the y = 0 plane. Thus, if the potential resulting 
from t his type of discontinuity (see fig. 32) is known, then 
the potential resulting from any combination of discontinu- 
ities of this type can be found. The potential for this type 
of distribution is denoted by </> 0 . By the use of cylindrical 
coordinates, as shown in figure 32, the potential at the point 
(x,p,d) can be expressed in terms of <£ n by 

<t>(x,p,d) = <l> 0 (x,p } d) — <t>o(x,p,0 — ^ ^ (64) 


z 



(O) 


z 



(b) 


z 


z 




Mach line from 
adjacent fin 


/--- Mach line from 
/ adjacent fin 



/ Mach line from 
root section 


Mach lines 


(e) 

(a) One fin. 
(c) Four fins. 


(b) Two fins. 

(d) Six fins. 

(e) Eight fins. 

Figure 29. — Regions of similar disturbances for tails consisting of 
rectangular fins. 

(50) and is given as follows: 

For ()<6<Z, 


<t>n(x,P ,0)= — j^xV-C 2 — 0V + 

-2 &P 


/3 2 p 2 ( 1 — 2 cos 2 0)ln 


z+V* 2 -0V 


o a a ■ -1 / /3p COS 6 \ 

2.rfip cos 6 sm ( - i— 

WP-0V sin 2 V 

20V cos # sin « tail" 1 (- ,)] 


Equation (64) follows from figures 31 and 32. 

The potential function <£ 0 was evaluated by use of equation 


(65a) 
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For — ^ 0 ^ 7r , 






0V+ 


0 2 p 2 ( 1 — 2 cos 2 0)ln 


Pp , 

X + V^ 2— i3 2 p 2 


The potential in region HI of the fin surface is a special 
case ^0=^J °f the preceding equation. Thus, the potential 
in region III is given (in Cartesian coordinates) by 


<t>(r, 0- tZ )=^. (x .sin 1 ; -fid n — -1^=^) 

\ ® x + ^x*—p z / 


=.) (67) 


_ . _ . , / Pp cos 0 \ . 

2 xpp cos 0 sin -1 ( o ■ ■ — )+ 

\^lx 2 -p 2 p 2 sin 2 0/ 

20 2 p 2 sin 0 cos 0 tan' 1 ( °* Y1 (6 

VV^ 2 -^ 2 p 2 /J 


5b) 


From equations (64) and (65), the potential function in the 
region of interaction, which is not affected by the tips, is 

given by the following equation for 


- 2 •*> , H . , *J_ r „.+ 

2xPp [~ cos 0 sin -1 ( Pp . cos ^ — V 

[_ \^x 2 —0 2 p 2 sin 2 0/ 


Sill 0 Sill 


yjx 2 —f3 2 p 2 sin 2 6j 

»>sm» \-| , 

\\x 2 — p 2 p 2 cos 2 0 / J 

_! jrjcot 0| 


20 2 p 2 cos 0 sin 0 ( tan 


tan 


_j a; | tan 0 

■yjx 1 


(* 

tan 0 | \) 

-(8V/J 


-\/x 2 —fi 2 p i 


( 66 ) 


/ 



Figure 30. — A cross section of the normal-velocity distribution on two 
plane surfaces representing the region of interaction for a rolling 
tail with four fins. 


From equation 
found to be 


(67) the pressure-difference 


n 8 pz . _! 02 
AC P =—f^— sin 1 — 
tVP x 


coefficient is 

( 68 ) 


z 



Figure 31. — A cross section of the velocity discontinuity distribution 
used to find the potential in part of the region of interaction for 
a tail of four fins. 


z 



Figure 32. — A cross section of the velocity discontinuity distribution 
associated with the function <£ 0 for a tail with four fins. 
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As previously stated, the potential in region IV (see fig. 29) 
is a combination of the potentials in regions I, II, and III. 
Assume that the fins are infinitely long. In this case only 
two regions (II and III) exist since regions I and IV are 
affected by the tip. The effect of the tip can be taken into 
account by adding a potential which has zero normal velocity 
on the fin and the negative of the pressure of the infinite fin 
in the plane of the fin outboard of the tip. The value of 
such a potential on the fin is given by the difference between 
the potential of region I and the potential of region II. (This 
potential is only the effect of the tip on a semi-infinite rolling 
wing.) Thus, the potential in region IV is the potential in 
region III plus the difference between the potential of region 
I and the potential of region II. Mathematically, the po- 
tential in region IV is given by 


<t>(x,0-,z)= =2- 


xz , pz . xz 

L~T C03 T+T 




KT _2_ 0V(i _z )( 2+ H) _2!ln i+v3 


(69) 


From equation (69) the pressure-difference coefficient is 
given by 



TAIL CONSISTING OF SIX FINS 

The pressure and potentials on the surface of the tail con- 
sisting of six fins can be obtained in a manner similar to that 
used for the tail consisting of four fins. The pressure and 
potentials in regions 1 and II (fig. 29 (d)) are the same as the 
pressure and potentials in the corresponding regions for tails 
consisting of one, two, or four fins. Regions III and IV are 
affected by the interaction between adjacent fins. The po- 
tential in region IV is made up of a combination of the po- 
tential in regions I, II, and 111; therefore, the main problem, 
as for the case of four fins, is the determination of the poten- 
tial in region III. 

The potential in region III is the same as the potential for 
a tail consisting of six infinitely long fins. The induced ve- 
locities normal to two of the planes of the fins are illustrated 
in figure 33. For two plane surfaces parallel to the stream 
direction and intersecting at an angle of 60°, the potential 
in region III can be obtained by a distribution of discontinui- 
ties in velocity as illustrated in figure 34. Note that the 
potential in region III can be made up of a combination of 
the potentials from a velocity discontinuity as shown in 
figure 35. The potential from this type of discontinuity is 
denoted by </>,,. By use of cylindrical coordinates as shown 
in figure 35, the potential at the point (x,p,d) can be expressed 
in terms of </> 0 by 

4>{r,p,6) = <t> n (x,p,6— ^)— <£„ ( x,p,6— )— <t>„ (/.p,0+^) (71) 


Equation (71) follows from figures 34 and 35. 

The potential function <£ 0 was evaluated by use of equation 
(50) and is given by the following equation for O^0^7r: 

<t> n {x,p,e)=- V ^~- 9 - {x-^p sin 9) (72) 

From equations (71) and (72) the potential function in the 
region of interaction, which is not affected by the tips, is 

given by the following equation for ^^0^^: 

2 

4>(r,p,0)=^y h 3(1 — 2 cos 2 0)+2 cos 0 sin 0] (73) 

The potential in region III of the fin surface is a special 
case ^ 0=7^) of equation (73). The potential in this region is 
given (in Cartesian coordinates) bv 

2 r 

<t>{x,Q-,z)= VZ ^ ] (74) 



Figure 33. — A cross section of the velocity distribution normal to the 
planes of two fins of a rolling tail consisting of six fins. 
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From equation (74) the pressure-difference coefficient in 
region III is found to be zero. 

The potential in region IV is a combination of the poten- 
tials in regions I, II, and III; this can be shown in the same 
way as the potential in region IV of the tail consisting of 
four fins was shown to be a combination of potentials from 
other regions. Specifically, the potential in region IV for 
the tail consisting of six fins is the potential in region 111 
plus the difference between the potentials of region I and of 
region II. Mathematically, the potential in region IV is 
given by 


0 ,z)= 


2 V 

TV 




1 

3 


(¥— )VG-*+H)J 


( 75 ) 


From equation (75) the pressure-difference coefficient is 
given bv 


A<7,= - 


Hp 

rp'V 




(76i 

TAIL CONSISTING OF EIGHT FINS 

The pressure and potential on the surface of the tail 
consisting of eight fins can be found by utilizing the potential 
functions </>„ used in finding the pressure and potentials on 
the surface of the tail consisting of four and six fins. The 
pressure and potentials in regions I and II are the same as 
the pressure and potentials for the corresponding regions of 
the other tails. The potentials in regions III, IV, V, and 
VI (see fig. 29(e)) are affected by the interaction between 
adjacent fins. Since the potentials in regions V and VI are 
combinations of the potentials in the remaining regions, the 
main problem is to find the potentials in regions III and IV. 

The potentials in regions III and IV are the same as the 
potential for a tail of eight infinitely long fins. The induced 
velocity normal to two of the planes of the fins is illustrated 
in figure 36. From the results for two plane surfaces inter- 
secting at an angle of 45°, the potentials in regions III and 
IV can be obtained by a distribution of discontinuities in 
velocity as illustrated in figure 37. The potential resulting 
from the distribution of discontinuities in velocity as illus- 
trated in figure 37 can be obtained from a distribution of 


z 



Figure 35. — A cross section of the velocity discontinuity distribution 
associated with the function <j> Q for a tail of six fins. 


! discontinuities as possessed by the potential function <£ 0 
used in connection with the four-finned tail. The potential 
function </> 0 used in connection with the four-finned tail 
was evaluated only in the region affected by the root 
sections of the fins. For the case of the eight-finned 
tail, interaction occurs between adjacent fins in regions 
which are not affected by the root sections of the fins. 

I The potential function </> 0 of the four-finned tail in the region 
; not affected by the root section must be known. In this 
region, the potential functions <£ 0 for the four- and six-finned 
tails are the same. 

From figures 32 and 37 the potential function in the region 
of fin interaction which is affected by the root sections of 

the fins is, for ^ 

<K* ,p,o) z =—<l>o(x,p,d)+<i>n (x,p,e—^\— 


<h (x,p,d— 0+4>o (a:,p,0+!) (77) 
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where is given by equations (65). 

The potential in the part of the region of interaction which is affected by the root section is (from eqs. (65) and (77)), 


for|g0g^ 
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(78) 


The potential in region III is a special case ^0=^ of equation (78). Setting 0==^ in equation (78) yields (in Cartesian 
coordinates) 

sin- S±-,j2sin- ^=J/ =jjp+U* tan \y£ r p) <”> 

From equation (79) the pressure-difference coefficient is found to be given by 


A (7 =.*££ /gin-1 £i_ 

p tt(3 c y x 


\ 2 sin -1 


02 


V 2x 2 — 


(80) 


By inspection of figures 35 and 37 the potential function in the region of fin interaction, which is not affected by the 


root sections of the fins, is, for ~ i£0 — 


4 i (x , p ,0) 0o (x, p .0) d - <f>n ^ x, p ,0 ^ ^ <f> 0 ( x, p ,0 ^ 0o (x, p ,0 ^ ) 


(81) 


where 0 O is given by equation (72). Substituting equation (72) into equation (81) yields (remember that 4>» is zero 

» > Jr 

upstream of the Mach cone from the y- axis), for — 0 5= y j 


4>(x,p 


,0)=- ~ ^(l — L) sin 0+' ° s ^J+pp 2 cos 0 sin 0 — 


(1 — 2 sin 2 0) — 


PP (cos 0 + si n0) 

dV 2 


(O: if a - < ^ (sin 0— cos 0) ) , n • , • „ 

) V 2 ' L^COS0 ° ; ^ PS1 

b— - p (sin0— cos0); ifxi> ^ P (sin0— cos0)( ^ (x-0P sin 0; if x>0p sin 0 

\ V 2 V 2 / 


(82) 


The potential in region IV is a special case ( 0=^) equation (82). Setting 9=— in equation (82) yields (in Cartesian 
coordinates) 

0(z,O-, 2 )=^[(1-V2)x+|8x] (83) 
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From equation (83) the pressure-difference coefficient is found to be given by 




(84) 


The potential in region V is the potential in region IV plus the difference between the potential in region I and the 
potential in region II; thus, from equations (62a), (62b), and (83), the potential in region V is found to be given by 




(85) 


From equation (85) the pressure-difference coefficient is found to be given by 

8p 
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7T|8 V 


. _,ir (2 *1, Ub \/ b,x\ 
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( 86 ) 


The potential in region VI is the potential in region III plus the difference between the potential in region I and the 
potential in region II; thus, from equations (62a), (62b), and (79), the potential in region VI is found to be given by 
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From equation (87), the pressure-difference coefficient ir. region VI is found to be given by 
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DISCUSSION OF RESULTS FOR ROLLING TAILS 

Illustrative plots of the chordwise and spanwise pressure distributions across one fin for tails with various numbers of 
fins are shown in figure 38. Figure 39 shows illustrative plots of the spanwise loadings on one fin for tails made up of 
various numbers of fins. 

The potential function <t> 0 used in finding the pressures and potentials for the tails consisting of four, six, and eight fins 
could be used in finding pressures and potentials for tails consisting of any even number of fins provided that the region of 
interaction between adjacent fins does not affect the tip. The restriction on the region of interaction causes the range of 
validity to decrease as the number of fins is increased. The range of validity could be extended, however, by use of a pressure 
potential cancellation method such as given in references 14 and 15. 

From the potential, the damping in roll per fin was calculated. Table I presents the results of these calculations. Figure 
40 presents the variation of the damping in roll per fin with A0 for tails made up of various numbers of fins. For a given 
Mach number (0 constant), figure 40 shows the variation of the damping with aspect ratio. Figure 41 presents the variation 
of the damping in roll per fin with Mach number for tails consisting of various numbers of fins with a fin aspect ratio of 1.5. 


Langley Aeronautical Laboratory, 
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Chordwise 


(a) Tail of one fin. 

Figure 38. — Chordwise and spanwise pressure distributions on a fin of aspect ratio 1.5 at a Mach number of \2 for tails 

consisting of rectangular fins. 
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(b) Tail of two fins. 
Figure 38. — Continued. 
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(c) Tail of four fins. 
Figure 38. — Continued. 
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Mach lines 



(d) Tail of six fins. 
Figure 38. — Continued. 
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(e) Tail of eight fins. 
Figure 38. — Concluded. 



Figure 39. — Spamvise loading on a fin of aspect ratio 1.5 at a Mach 
number of V 2 for tails consisting of rectangular fins. 
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Mach number 


Figure 41. — An illustrative variation of Ci p per fin with Mach number for tails consisting of rectangular fins with a fin aspect ratio of 1.5. 
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TABLE T 


DAMPING-IN-ROLL COEFFICIENT PER FIN 


Number qq per fin Valid for 

of fins j ” 

1 

1 + 4.4/3-24.1*02 + 32/l 3 0 3 
24 

/10S1 

2 

1 +8/10-48/l 2 0 3 + 64/l 3 0 3 
48 

Aft *i 

4 

" J 0 3 (1 + M0— 48^0*+64/l 3 0 3 )] 

/ 10^1 

6 

-.i' [-,3s + ns a 

All 

<0. 

1 8 

-.lU-iH+i bO+w-w+iww] 

/10S/2 
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